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Finite  Deformations  and  Internal 
Forces  in  Elastic-Plastic  Crystals: 
Interpretations  From  Nonlinear 
Elasticity  and  Anharmonic  Lattice 
Statics 

Large  deformation  kinematics  and  internal  forces  arising  from  defects  in  crystalline 
solids  are  addressed  by  a  nonlinear  kinematic  description  and  multiscale  averaging 
concepts.  An  element  of  crystalline  material  with  spatially  uniform  properties  and  con¬ 
taining  defects  such  as  dislocation  lines  and  loops  is  considered.  The  average  deforma¬ 
tion  gradient  for  this  element  is  decomposed  multiplicatively  into  terms  accounting  for 
effects  of  dislocation  flux,  recoverable  elastic  stretch  and  rotation,  and  residual  elastic 
deformation  associated  with  self-equilibrating  internal  forces  induced  by  defects.  Two 
methods  are  considered  for  quantifying  average  residual  elastic  deformation:  continuum 
elasticity  and  discrete  lattice  statics.  Average  residual  elastic  strains  and  corresponding 
average  residual  elastic  volume  changes  are  negligible  in  the  context  of  linear  elasticity 
or  harmonic  force  potentials  but  are  not  necessarily  inconsequential  in  the  more  general 
case  of  nonlinear  elasticity  or  anharmonic  interactions.  [DOI:  10.1115/1.3183773] 

Keywords:  nonlinear  elasticity,  plasticity,  lattice  statics,  dislocations,  multiscale 
modeling 


1  Introduction 

By  typical  definition,  plastic  deformation  in  crystalline  solids 
takes  place  by  the  motion  or  flux  of  distributions  of  dislocations 
[1,2].  An  element  of  material  through  which  a  net  flux  of  disloca¬ 
tions  has  passed  exhibits  a  permanent  (plastic)  shape  change  [3]. 
The  dislocation  lines  within  the  element,  i.e.,  displacement  dis¬ 
continuities  across  the  slip  plane  in  the  context  of  Volterra  defects 
in  elastic  continua,  may  also  contribute  to  this  plastic  deformation 
[4-7],  However,  plastic  deformation  in  this  context  does  not  ex¬ 
plicitly  account  for  the  additional  change  in  dimensions  of  the 
body  due  to  residual  elastic  deformation  associated  with  local 
stress  fields  (i.e.,  eigenstresses  [7])  induced  by  defects.  For  ex¬ 
ample,  dislocation  glide  preserves  the  volume  of  the  crystal,  and 
tangential  displacement  (i.e.,  slip)  discontinuities  do  not  alter  the 
volume  occupied  by  the  material.  Yet  ample  evidence  suggests 
that  dislocation  lines  affect  the  volume  of  crystals  [8-11].  Stored 
energies  associated  with  elastic  fields  of  defects  are  important 
because  they  affect  recrystallization  [9,12]  and  the  fraction  of 
stress  power  converted  to  temperature  rise  at  high  deformation 
rates  [13]  that  can  lead  to  strain  softening  and  shear  localization  in 
metals  during  dynamic  failure  events.  Large  numbers  of  disloca¬ 
tions,  twins,  and  stacking  faults  can  be  generated  during  shock 
loading  [14],  and  presumably  the  corresponding  volume  changes, 
shape  changes,  and  stored  energies  associated  with  local  stress 
fields  of  these  defects  affect  the  observed  response  (e.g.,  pressure- 
volume  or  stress-strain  profiles)  under  such  conditions  [15]. 

Defects  of  interest  in  the  present  work  are  those  that  can  be 
described  in  a  continuum  sense  by  closed  displacement  disconti¬ 
nuities  tangential  to  an  internal  surface  (i.e.,  crystallographic 
plane)  in  a  volume  element  of  material.  These  include  gliding 
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straight  or  curved  dislocations  lines,  dislocation  loops,  as  well  as 
partial  dislocations.  Motion  of  dislocation  defects  through  a  re¬ 
gion  of  the  crystal  results  in  plastic  shape  change  as  mentioned 
above  but  preserves  the  lattice  spacing  [3];  this  requires  coopera¬ 
tive  motion  of  leading  and  trailing  partials  for  the  case  of  partial 
dislocations.  Residual  elastic  stress  fields  from  disclination  lines 
and  loops  are  also  considered  [6,16,17],  though  their  contribution 
to  the  plastic  shape  change  [18]  is  not  addressed  explicitly.  Also 
not  considered  are  voids,  fractures,  and  point  defects  [19,20]  that 
may  have  finite  volume  and  affect  the  lattice  differently.  Disloca¬ 
tion  climb  is  also  not  considered  here,  because  it  generally  in¬ 
volves  vacancy  migration  [21], 

Modeling  efforts  directed  toward  large  deformation  kinematics 
of  crystalline  solids  feature  a  long  history.  In  the  absence  of  de¬ 
fects,  thermal  or  electromagnetic  effects,  or  internal  relative  dis¬ 
placements  among  atoms  of  different  sublattices,  a  deformation 
gradient  F  applied  uniformly  over  an  element  of  mechanically 
stressed  material  and  the  linear  mapping  FE  of  its  interatomic 
bond  vectors  (i.e.,  primitive  Bravais  lattice  vectors  and  basis  vec¬ 
tors  at  each  lattice  point)  coincide.  This  is  typically  called  the 
Cauchy-Bom  hypothesis  [22,23]  and  is  written  as  F=FE.  When 
plastic  deformation  takes  place,  the  description  is  typically  ex¬ 
tended  to  F  =  FeFp  [3,24,25],  where  Fp  accounts  for  the  presence 
and  motion  of  dislocations  or  other  defects.  The  idea  of  a  relaxed 
intermediate  or  natural  configuration  implicit  in  this  multiplicative 
decomposition  was  introduced  somewhat  earlier  [26].  Early  stud¬ 
ies  of  continuously  distributed  lattice  defects  considered  self¬ 
stresses  and  residual  elastic  deformations  resulting  from  defects 
[27-29].  In  such  instances  the  plastic  deformation  Fp  inducing 
discontinuities  in  the  first  place  was  often  not  considered  explic¬ 
itly,  nor  was  elastic  deformation  FE  associated  with  external  load¬ 
ing,  since  the  element  of  material  containing  defects  was  in  a  state 
of  self-stress.  Defect  content  associated  with  heterogeneous  elastic 
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and  inelastic  deformations  (i.e.,  tangent  mappings  in  the  multipli¬ 
cative  decomposition)  can  be  described  using  tensor  calculus  and 
differential  geometry  [3,16,24,30-37], 

In  this  work  a  three-term  decomposition,  F=FEFIFP,  is  as¬ 
signed  to  describe  the  kinematics  of  a  volume  element  of  crystal¬ 
line  material  containing  lattice  defects.  For  simplicity,  isothermal 
conditions  are  assumed  (e.g.,  no  thermal  expansion).  Elastic  tan¬ 
gent  map  Fe  used  here  is  analogous  to  that  of  Bom,  Huang,  and 
Ericksen  [22,23,38],  Plastic  tangent  map  Fp  results  from  cumula¬ 
tive  effects  of  fluxes  of  dislocations  [1,2,39]  and  slip  discontinui¬ 
ties  associated  with  defects  within  the  element  [4,5],  Intermediate 
term  F1  (or  its  stretch  or  rotational  components)  has  been  intro¬ 
duced  in  a  number  of  works  with  various  definitions 
[16,34-36,40-42],  Here  F1  represents  the  average  residual  elastic 
deformation  of  the  element  induced  by  local  stress  and  strain 
fields  of  defects  contained  within,  and  corresponds  to  the  volume 
average  of  residual  elastic  deformation  in  the  context  of  defect 
field  theories  [27-29], 

In  Sec.  3,  the  externally  unloaded  volume  element  is  treated  as 
an  elastic  body  with  homogeneous  elastic  constants,  in  static  equi¬ 
librium,  and  containing  internal  displacement  discontinuities 
across  which  traction  is  continuous.  The  element  is  free  of  exter¬ 
nal  traction  but  may  support  residual  stresses.  The  self¬ 
equilibrium  conditions  result  in  an  integral  equation  for  F1.  For 
linear-elastic  constitutive  behavior,  F1  reduces  to  the  unit  tensor  in 
rectangular  Cartesian  coordinates,  but  F1  does  not  so  reduce  for 
nonlinear  behavior.  Following  previous  studies  [29,43,44],  vol¬ 
ume  changes  resulting  from  the  stored  energy  of  defect  lines  are 
considered  for  cubic  and  isotropic  crystals.  A  proportionality  rela¬ 
tionship  between  the  line  density  of  dislocations  or  disclinations 
per  unit  reference  volume  of  the  solid  and  this  residual  volume 
change  is  obtained.  Analytical  elasticity  solutions  for  defect  ener¬ 
gies  are  considered  along  with  experimental  data  [45]  to  confirm 
accuracy  of  the  theory. 

In  Sec.  4,  the  externally  unloaded  volume  element  is  identified 
as  a  set  of  discrete  atoms  in  static  equilibrium  and  free  of  external 
forces.  According  to  this  description,  all  interatomic  forces  vanish 
in  a  perfect  lattice  free  of  defects  (e.g.,  no  bond  stretching),  but 
nonzero  forces  may  exist  among  atoms  when  they  are  arranged  in 
an  imperfect  way  (i.e.,  when  defects  are  present).  In  this  case,  the 
self-equilibrium  conditions  are  equivalent  to  vanishing  of  static 
components  of  the  average  virial  stress  for  the  set  of  atoms 
[38,46],  The  potential  energy  is  expanded  about  a  perfect  refer¬ 
ence  state,  leading  to  expressions  for  stress  and  elastic  stiffness 
tensors  in  terms  of  harmonic  and  anharmonic  parts  of  the  poten¬ 
tial.  When  atomic  interactions  are  harmonic,  equilibrium  demands 
that  F1  as  defined  here  should  reduce  to  the  unit  tensor,  analo¬ 
gously  to  the  linear-elastic  continuum  interpretation.  However,  in 
the  more  general  scenario  of  anharmonic  interactions  [47-49]  that 
may  be  physically  significant  in  the  vicinity  of  defect  cores  where 
large  atomic  displacements  arise,  the  derivations  suggest  F1  could 
be  non-negligible. 

A  number  of  theoretical  models  have  introduced  a  term  akin  to 
F1  to  represent  elastic  fields  of  crystal  defects  [16,34,35,41], 
though  previous  works  did  not  always  include  precise  mathemati¬ 
cal  definitions  for  F1  or  quantitative  estimates  of  its  magnitude  in 
engineering  materials.  The  present  work  gives  exact  definitions 
for  F1  using  anisotropic  nonlinear  continuum  elasticity  (Sec.  3) 
and  discrete  lattice  statics  (Sec.  4).  Also,  the  present  work  dem¬ 
onstrates  how  residual  elastic  volume  changes  (~det  F1)  can  be 
estimated  for  defect  lines,  extending  earlier  analysis  [43]  to  ad¬ 
dress  geometric  nonlinearity.  Previous  authors  [8,10,43]  did  not 
emphasize  possible  shape  changes  associated  with  F1  (though 
Toupin  and  Rivlin  [44]  mentioned  shape  changes  briefly)  nor  did 
they  attempt  to  place  their  derivations  in  the  context  of  multipli¬ 
cative  elastoplasticity  as  demonstrated  in  the  present  work.  The 
new  atomic-scale  definition  of  F1  in  Sec.  4  may  prove  useful  in 
the  context  of  multiscale  computations  of  nonlinear  elastic  and 
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Fig.  1  Configurations,  surface  coordinates,  traction  vectors, 
and  tangent  mappings 


plastic  properties  of  crystals  [50,51]  and  offers  a  more  accurate 
treatment  of  defect  cores  than  classical  continuum  elasticity. 

The  following  notation  is  used.  Scalars  and  individual  compo¬ 
nents  of  vectors  and  tensors  are  written  in  italics,  while  vectors 
and  tensors  are  written  in  bold.  Einstein’s  summation  convention 
applies  for  repeated  indices.  The  •  symbol  denotes  the  scalar  prod¬ 
uct  of  vectors  (a'b=aaba=alb\  +a1b2+aib^),  while  ®  indicates 
the  outer  product  ((a<8>b )ab=aabb).  Juxtaposition  of  second-rank 
tensors  implies  summation  over  one  set  of  adjacent  indices 
((AB )ac=AabBbc).  Indices  in  parentheses  are  symmetric  (2 
=Aab+Aba).  Superposed  —  1  denotes  inversion.  Subscripted  com¬ 
mas  denote  partial  differentiation. 

2  Kinematics  of  Deformation 

Consider  a  volume  element  of  a  crystal  containing  lattice  de¬ 
fects.  This  element  is  of  a  size  larger  than  the  atomic  spacing  but 
may  be  of  smaller  dimensions  than  the  entire  crystal.  Static  and 
isothermal  conditions  are  assumed.  Deformation  of  the  element  is 
described  by  mappings  between  tangent  spaces  of  configurations, 
as  shown  in  Fig.  1 .  Reference  configuration  B0  is  a  perfect  lattice. 
Position  vector  of  atom  m  (e.g.,  the  equilibrium  position  of 
its  nucleus)  relative  to  the  origin  of  a  Cartesian  coordinate  system 
in  an  infinite  lattice  in  the  reference  configuration  is  given  by  the 
periodicity  relations  [22,23] 

3 

R<m>  =  z{m)  +  2  i(m)k^k  ( 1 ) 

k=\ 

where  i(m)b  are  integers  and  A*  are  the  primitive  Bravais  lattice 
vectors.  Basis  vector  z denotes  the  reference  position  of  atom  m 
of  the  basis  at  each  point  on  the  Bravais  lattice  relative  to  the 
origin  of  the  corresponding  primitive  unit  cell.  For  a  monatomic 
basis,  Z(m)  =  0. 

2.1  Plastic  Deformation.  Intermediate  configuration  B  in  Fig. 
1  differs  from  B0  due  to  influences  of  cumulative  motion  of  lattice 
defects  and  perturbations  of  atomic  positions  resulting  from  these 
defects.  Configuration  B  is  by  definition  free  of  external  traction 
(F=0)  and  free  of  internal  stresses.  For  a  crystal  containing  a 
single  dislocation,  a  singular  plastic  deformation  map  can  be  de¬ 
fined  as  [4,5] 

Fp  =  l  +  b<8>M<5(X*M)x  (2) 

where  1  is  the  unit  tensor,  b  is  the  Burgers  vector,  M  is  the  normal 
to  the  slip  surface  in  the  reference  configuration,  8  is  the  Dirac 
delta  function,  and  \  is  the  characteristic  function  that  is  unity  at 
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reference  coordinates  X  on  slipped  surface  X  and  zero  elsewhere. 
The  Cartesian  coordinate  system  for  X  is  chosen  with  its  origin  in 
the  slip  plane.  A  dislocation  density  tensor  corresponding  to  Eq. 
(2)  is  [4] 

a  =  b®g0S(L)  (3) 

where  g0  is  the  unit  tangent  to  the  dislocation  line  L  in  the  refer¬ 
ence  configuration.  From  Eqs.  (2)  and  (3), 


Relation  (9)  applies  for  straight  dislocations;  similar  expressions 
can  be  formulated  for  dislocation  loops  [4],  Introducing  the  nota¬ 
tion  g1  Xv'  =  m'v',  u'  =  |v1|,  b'  =  b's\  and  ft'  =  |b'|,  where  m'  is  the 
unit  normal  to  the  slip  plane  and  s’  is  a  unit  vector  in  the  direction 
of  slip,  Eq.  (9)  reduces  to  [5,53] 

LP  =  S  p‘b  W  ®  m'  =  2  y's1'  ®  m'  (10) 


b  =  -  I  FprfX=  aNdA  (4) 

J  x*  Ja 

where  the  line  integral  in  the  first  of  Eq.  (4)  takes  place  across 
coordinates  above  (X+)  and  below  (X-)  the  slipped  surface,  and 
integration  in  the  second  of  Eq.  (4)  proceeds  over  oriented  area  A 
with  unit  normal  N  such  that  g0  •  N  =  1  at  the  intersection  of  L  and 
A.  If  defects  are  introduced  sequentially,  a  logical  nonsingular 
extension  of  Eq.  (2)  for  a  volume  element  containing  multiple 
dislocation  lines  is 


Fp=n 

m=  1 


b  (8)  M<5(X)rfF  =  11  (1  +  V“’b  ®  MX) 

/  m=  1 


Rate  equation  y  =  p'b'vl  [1]  implies  that  only  mobile  dislocations 

contribute  to  Eq.  (10).  Since  s' ±m‘,  jp=Jp  tr  Lp=0  and  Jp=l.  In 
summary,  according  to  definitions  in  the  context  of  single  dislo¬ 
cations  (2)  and  (5),  dislocation  flux  (9),  or  crystal  plasticity  theory 
(10),  plastic  deformation  is  isochoric. 

2.2  Lattice  Deformations.  Let  FE  denote  the  tangent  map  to 
externally  stressed  spatial  configuration  B  (generally  with  t  A  0  in 
Fig.  1)  from  self-equilibrated  configuration  B  (t=0  in  Fig.  1).  Let 
a*  denote  the  primitive  Bravais  lattice  vectors  of  Eq.  (1)  mapped 
to  B : 


a‘=  FTA* 


(ID 


(5) 

where  the  product  is  taken  over  n  dislocations,  each  with  a  possi¬ 
bly  different  Burgers  vector  b  and  constant  normal  vector  M  to  an 
assumed  flat  slipped  surface  of  reference  area  X.  The  volume  of 
the  element  in  the  reference  configuration  is  denoted  by  V.  Note 
that  Eq.  (5)  depends  on  the  order  in  which  each  dislocation  is 
introduced.  Since  blM,  Eq.  (5)  consists  of  a  product  of  simple 
shears  and  volume  is  conserved,  i.e.,  7p  =  det  Fp=l.  For  a  collec¬ 
tion  of  j  straight  stationary  dislocation  line  populations  of  density 
p'=  lim(L'/ V),  Burgers  vector  b',  unit  tangent  g'0,  and  plane  nor- 

mal  M',  the  dislocation  density  in  Eq.  (3)  and  Burgers  integral  in 
Eq.  (4)  can  be  extended  as 

a=\T1  f  b  ®  g0dL  =  2  p'b'  ®  gf)  (6) 

J  L  1=1 


B  =  aJ\dA  =  2  p'b'(£o  *  NcM)  (7) 

J  A  J  A  1=1 

with  B  the  summed  projection  of  all  local  Burgers  vectors  b  of 
dislocation  lines  passing  through  area  A.  The  line  integration  in 
Eq.  (6)  proceeds  over  all  dislocations  of  total  length  L.  The  second 
of  Eq.  (6)  pertains  to  j  dislocation  populations  with  the  same 
tangent  line  and  local  Burgers  vector  for  each  value  of  i. 

A  more  general  definition  of  average  plastic  deformation  of  a 
volume  element  accounts  for  the  history  of  generation  and  motion 
of  defects  within  it,  in  which  case  Fp  is  defined  as  the  solution  of 

Fp  =  LpFp,  FPU  =  1  (8) 

Average  plastic  velocity  gradient  Lp  is  dictated  by  the  flux  s  of 
mobile  dislocations  [2,17,39,52]: 

LPo,  =  z*xSgxSfh  ?*x«  =V~I  f  bir‘g‘xv‘8dL  =  2  pV^W8 

Jl  1=1 

(9) 

where  e  is  the  permutation  tensor,  and  g'  and  v1  are,  respectively, 
the  uniform  tangent  line  and  velocity  of  every  dislocation  in  popu¬ 
lation  i.  All  vector  and  tensor  indices  written  in  Greek  font  in  Eq. 
(9)  are  subject  to  the  Einstein  summation  convention  and  are  re¬ 
ferred  to  configuration  B.  Since  dislocation  segments  have  perpen¬ 
dicular  velocities  and  tangent  lines,  JP=JP  trLp=0  and  Jp=  1. 


with  1  the  shifter  between  coordinate  frames  in  B0  and  B.  When 
coincident  coordinate  systems  are  used  in  these  two  configura¬ 
tions,  1  =  1.  Note  that  lattice  vectors  are  not  affected  by  Fp,  in 
accordance  with  continuous  dislocation  theory  [3,34]  and  crystal 

plasticity  theory  [53].  Lattice  vectors  deform  from  B  to  B  accord¬ 
ing  to 


a*  =  FEa*  =  FLlA*  (12) 

where  the  total  lattice  deformation  for  the  volume  element  is  FL 
=  FeF*.  Under  a  homogeneous  deformation,  basis  vectors  z ^  in 
Eq.  (1)  deform  similarly  to  the  primitive  Bravais  lattice  vectors  A* 
in  Eq.  (12),  though  the  final  position  of  atom  m  may  also  change 
by  an  additional  translation  relative  to  its  neighbors  [22],  This 
translation  may  be  due,  for  example,  to  polarization  of  a  dielectric 
crystal  in  an  applied  electric  field  or  inner  displacements  among 
sublattices  in  a  noncentrosymmetric  polyatomic  crystal  [54],  Re¬ 
lations  (11)  and  (12)  generalize  the  Cauchy-Born  hypothesis 
[22,23],  distinguishing  the  effects  of  recoverable  (FE)  and  residual 
(F1)  deformations  on  the  lattice,  and  reduce  to  the  prescription  of 

Ericksen  [23]  when  the  effects  of  defects  are  absent  in  B,  such  that 
F-l.  Rigid  body  rotations  of  the  element  are  included  in  the 
rotational  part  of  the  polar  decomposition  of  FL,  denoted  by  RL. 
Reference  Bravais  lattice  vectors  A*  are  uniform  in  a  perfect  ho¬ 
mogeneous  lattice.  However,  when  defects  are  contained  within 
the  volume,  deformed  primitive  lattice  vectors  aA  and  a*  represent 
suitable  averages  of  local  lattice  vectors  a'A  and  a'k  that  may  not 
be  spatially  constant: 


Configuration  B  can  be  obtained  from  the  spatial  configuration  by 
cutting  the  volume  element  out  of  the  stressed  body  thereby  re¬ 
lieving  the  external  traction,  relaxing  any  possible  internal  viscous 
and  inertial  forces,  and  then  rotating  this  element  by  RE_1. 

2.3  Total  Deformation.  Total  deformation  gradient  F  for  the 
volume  element  is  defined  by  the  surface  integral  [55,56] 

F  =  V^1  J  x®N dS  (14) 

J  s 

where  x  are  spatial  coordinates  of  the  deformed  image  of  refer¬ 
ence  surface  S  enclosing  the  volume  element  with  unit  outward 
reference  normal  N,  as  shown  in  Fig.  1.  The  element  may  contain 
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discontinuities  in  internal  displacement  and  gradients  of  displace¬ 
ment;  when  discontinuities  are  absent  and  after  Gauss’s  theorem  is 
applied,  Eq.  (14)  reduces  to  F°A=  \T1fxa^dV,  where  xa=xa(XA  ,t) 
is  now  smooth  within  V.  Assume  that  F  is  imposed  on  an  element 
of  material  via  Eq.  (14),  and  that  decomposition 

F  =  FeFiFp  =  FlFp  (15) 

applies,  as  implied  by  Fig.  1.  Plastic  deformation  Fp  is  known 
from  integration  of  Eq.  (8)  with  Eq.  (9)  or  Eq.  (10),  presuming 
kinetic  laws  are  available  for  the  dislocation  flux  in  Eq.  (9)  or  the 
slip  rates  on  each  slip  system  in  Eq.  (10).  Under  isothermal  con¬ 
ditions,  the  crystal  responds  to  applied  loading  elastically  such 
that 

2  =  2(UE)^>UE  =  UE(2)  (16) 

where  stretch  UE  =  RE_1FE  =  (FE7FE)1/2  is  an  invertible  function 
of  conjugate  stress  measure  2  that  vanishes  when  traction  t  =  0 
along  deformed  surface  S  of  the  element  in  Fig.  1.  The  particular 

form  of  Eq.  (16)  will  depend  on  the  orientation  of  the  lattice  in  B 
because  of  anisotropy.  Henceforth,  it  is  assumed  that 

Fl  =  ReUeF1  =  RlUeU\  F1  =  U1  (17) 

meaning  F1  is  a  stretch  (symmetric  with  six  independent  entries  in 
covariant  Cartesian  coordinates)  and  all  lattice  rotation  is  embed¬ 
ded  in  Re.  The  objective  of  Secs.  3  and  4  of  the  present  paper  is 
suggestion  of  approaches  to  obtain  F1.  Then  if  Fp  is  also  known  at 
a  particular  instant,  FE=FFP_1FI_1  can  be  found  from  Eq.  (15) 
and  the  average  external  stress  supported  by  the  element  can  be 
updated  according  to  constitutive  relation  (16). 


3  Nonlinear  Elastic  Approach 

The  dimensional  changes  of  a  nonlinear  elastic  body  in  a  state 
of  self-stress,  i.e.,  a  self-equilibrated  body  with  internal  residual 
stresses  but  no  traction  applied  to  its  external  boundaries,  are  de¬ 
rived  in  Secs.  3. 1-3.3.  The  body  may  contain  one  or  more  internal 
surfaces  across  which  traction  is  continuous  but  tangential  dis¬ 
placements  are  not.  The  treatment  is  specialized  in  Sec.  3.4  to 
address  volume  changes  in  cubic  crystals  and  then  isotropic  ma¬ 
terials.  Formulas  for  volume  changes  attributed  to  local  elastic 
stress  fields  of  line  defects  are  derived  in  the  isotropic  approxima¬ 
tion  in  Sec.  3.5.  Derivations  and  discussion  in  Secs.  3. 1-3.5  con¬ 
solidate  and  extend  prior  work  of  a  number  of  authors 
[8,29,43,44],  Additional  analysis  and  examples  follow  in  Secs.  3.6 
and  3.7. 

3.1  Average  Stress  Measures.  The  local  balance  of  linear 
momentum  is  written  in  rectangular  Cartesian  coordinates  as 

P“aa  +  B“  =  P{A“  (18) 

with  Ba  the  body  force  per  unit  reference  volume,  p0  the  reference 
mass  density,  and  Aa  the  material  acceleration.  In  reference  coor¬ 
dinates  XA,  Eq.  (18)  provides  the  relation 

( XAPaB)B  =  P“A  +  (PoAa-Ba)XA  (19) 

Consider  a  body  of  reference  volume  V  with  external  surface  S. 
The  body  may  contain  closed  internal  surfaces  across  which  dis¬ 
placement  from  the  reference  state  (e.g.,  a  perfect  lattice)  and 
stress  fields  are  discontinuous,  but  traction  per  unit  reference  area 
fg  across  internal  surfaces  is  continuous: 

toa  -  fQa  =  (P+“A  -  P~aA)MA  =  0  (on  2)  (20) 

where  +  and  —  denote  the  limiting  values  of  a  quantity  near  the 
surface  as  the  surface  is  approached  from  either  side,  MA  is  a 
normal  vector  to  an  internal  surface,  and  2  denotes  the  union  of 
such  internal  surfaces.  The  source  of  the  displacement  disconti¬ 
nuities  across  2  is  arbitrary  in  Eq.  (20);  however,  for  the  particular 
case  of  dislocation(s)  within  V,  jump(s)  in  displacement  across  the 


slip  planes  comprising  2  are  attributed  to  the  Burgers  vector(s) 
introduced  in  Eqs.  (2)  and  (5).  Integrating  Eq.  (19)  over  V  and 
applying  the  divergence  theorem. 


PaAdV=  XAP“BNBdS  +  XA(P+aB  -  P~“B)MBdt 

J  v  J  s  J  s 


+  Xa(B“  -  PoA“)dV 

J  v 


(21) 


where  NB  are  components  of  the  external  normal  to  V  of  Eq.  (14). 
Applying  Eq.  (20),  and  considering  now  a  body  in  static  equilib¬ 
rium, 


P"AdV  = 

J  v  J  s 


XAP 


" BNBdS  = 


(22) 


For  a  self-equilibrated  body,  fg=0  by  definition,  and  Eq.  (22)  re¬ 
duces  to 


V1  P“AdV  =  0  (23) 

J  v 

meaning  that  the  integrated  or  volume-averaged  first  Piola- 
Kirchhoff  stress  vanishes  over  the  reference  volume.  Now  con¬ 
sider  the  balance  of  linear  momentum  in  the  spatial  configuration, 

aabb  +  b“  =  Paa  (24) 

with  oab  =  dst(XAa)xaAPbA,  ba  the  body  force  per  unit  current  vol¬ 
ume,  p=p0  det(XAn),  and  aa  the  spatial  acceleration.  Spatial  ana¬ 
logs  of  Eqs.  (19)— (21)  are 

(xVc),c  =  <fh  +  (pa“  -  ba)xh  (25) 

t+a  -  t~a  =  (<T+ah  -  <T~ah)mh  =  0  (on  a)  (26) 


|  r i°bdv  =  j  xba“cncds  +  j  xh(cr+ac  -  cr  ac)mcdcr 

J  v  J  s  a 

+  xh(ba  -  Paa)dv 

J  II 


(27) 


with  x“  the  spatial  coordinates,  t“  the  traction  per  unit  current  area, 
< t  the  union  of  closed  internal  boundaries,  ma  the  unit  normal  to 
internal  surfaces,  na  the  unit  normal  to  external  surface  S  (Fig.  1), 
and  u  the  current  volume  of  the  body  enclosed  by  s.  For  a  body  in 
static  equilibrium. 


J  i/'hdv  =  J  xh(f,cncds  =  J  fxbds 

J  V  J  S  J  s 


and  for  a  self-equilibrated  body, 


J  V 


a“bdv  =  0 


(28) 


(29) 


i.e.,  the  average  Cauchy  stress  vanishes  over  the  spatial  volume  of 
the  body. 

3.2  Hyperelasticity.  In  the  present  treatment,  the  constitutive 
response  of  the  material  is  assumed  hyperelastic,  with  effects  of 
temperature  change  neglected.  Let  'F0=1^o(^'Ab)  denote  the  strain 
energy  per  unit  reference  volume  of  the  solid,  with  the  symmetric 
Lagrangian  strain 
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^ AB  ~  2 .,A°abX -,B  ~  °AB )  ~  °.A 


Au(a,B )  +  2  Uc,AU.,B 


(30) 


where  coincident  Cartesian  coordinates  are  assigned  in  reference 
and  spatial  configurations  and  ua=xa-  8‘AXA  is  the  local  displace¬ 
ment.  Let 


%  =  ~ CabcdEabEcd  +  ±-CabcdefE, 


abEcdEef+  (31) 


where  second-  and  third-order  elastic  constants  are  evaluated  from 
a  series  expansion  of  the  energy  density  about  the  unstrained  state 
of  the  crystal  as,  respectively, 


C 


■ABCD  _ 


^0 


c 


-'ABCDEF  _ 


<^0 


E=0 

(32) 

Usual  symmetries  of  elastic  moduli  are  evident  from  Eqs.  (31)  and 
(32).  The  local  first  Piola-Kirchhoff  stress  following  from  Eq. 
(31)  is 


dXaA 


du. 


a^A 


3dEAR 


\ ,B 


CA 


’+2C 


r  r  17 


(33) 


Using  symmetry  properties  of  the  moduli  from  Eq.  (32)  and  ne¬ 
glecting  terms  of  degree  higher  than  3  in  the  displacement  gradi¬ 
ents  uAB=  ^AuBB,  Eq.  (31)  can  be  written  as  [44] 


M/  _  lpABCDu  IpABCD  E 

Tu-  UA,BUC,D  +  2^  UA,BUE,CU.,D 


+  gC 


LrABCDEF , 


UA,BUC,DUE,F 


(34) 


3.3  Average  Residual  Elastic  Deformation.  Consider  exter¬ 
nally  unloaded  configuration  of  the  volume  element  labeled  B  in 
Fig.  1.  This  configuration  corresponds  to  the  deformed  but  self- 
equilibrated  body  obeying  Eqs.  (23)  and  (29).  The  reference  con¬ 
figuration  with  internal  displacement  discontinuities  described  in 
Eqs.  (18) — (34)  now  corresponds  to  stress-free  intermediate  con¬ 
figuration  B  of  Fig.  1.  As  discussed  by  De  Wit  [6],  displacement 
discontinuities  attributed  to  Volterra  dislocations  lead  to  disconti¬ 
nuities  in  local  plastic  strain  and  rotation  fields  across  slip  planes 
or  internal  surfaces  2  within  the  volume,  which  owing  to  the 
isochoric  nature  of  slip  is  preserved  according  to  V=JPV=V.  Such 
effects  are  quantified  according  to  Eqs.  (2)-(10),  wherein  discon¬ 
tinuities  are  introduced  by  the  plastic  deformation  Fp  between 
configurations  BQ  and  B.  Residual  deformations  are  introduced 
within  the  element  from  the  mapping  between  configurations  B 
and  B.  Locally,  such  deformations  are  treated  as  continuous,  dif¬ 
ferentiable,  and  elastic  between  B  and  B ,  with  corresponding  local 
displacement  gradient  and  Lagrangian  strain  fields  (i.e.,  eigen- 
strains)  within  the  element  denoted  by  uA  B  and  Eab>  respectively. 
This  implies  that  the  body  has  been  resealed  after  the  introduction 
and  passage  of  defects  in  agreement  with  the  definition  of  a  Vol¬ 
terra  or  Somigliana  dislocation  [18],  and  residual  elastic  deforma¬ 
tions  manifest  within  a  body  that  can  be  treated  as  locally  con¬ 
tinuous  in  B.  The  corresponding  definition  for  the  average  residual 
elastic  deformation  gradient  in  Eq.  (17)  is,  in  Cartesian  coordi¬ 
nates, 


F1  - 

rAB  - 


i 


X(A,B)‘ 


4V=  SAR  +  V 


1 


l(A,B)‘ 


dV 


(35) 


Relation  (35)  is  a  key  definition  that  will  be  used  repeatedly  in 
what  follows.  Mapping  the  self-equilibrium  conditions  in  Eq.  (29) 
to  the  reference  configuration  and  substituting  from  Eq.  (33), 


Z'A  f  <fbdu  =  #a#b\  xbcP“cdV=  [  (# 

J  „  J  V  J  V 


0 


■c  +  u  c)~ -dy=0 


du 


A,C 


(36) 

The  integral  in  Eq.  (36)  becomes,  upon  substitution  of  Eq.  (34) 
with  omission  of  terms  of  orders  higher  than  3, 

0=  f  (SBc  +  uBc)^dV  =  (  CABCDuCDdV+  f  (Caecdube 
J  V  ’  duA,C  J  v  J  v 

iJ. 


+  Cbecduae)uc  Ddv  +  ^  I  CABCDuECuhDdV 


+  ~  I  C ABCDEFuc  DuEFdV 


(37) 


Neglecting  products  of  order  2  in  displacement  gradients  in  a 
linear-elastic  body  with  spatially  constant  moduli,  Eq.  (37)  re¬ 
duces  to 


CABCD  uCMdV=  0^  u(CD)dV=0  —*  Fed  =  SCD  (38) 

J  v  -J  V 

since  CABCD  is  assumed  positive  definite.  Linear  elastic  approxi¬ 
mation  (38)  states  that  no  average  elastic  shape  or  volume  change 
occurs  in  a  homogeneous  self-equilibrated  body.  The  latter  would 
be  true  even  if  local  elastic  dilatation  from  defects  does  not  van¬ 
ish,  as  is  exhibited  in  straight  line  edge  dislocation  and  wedge  and 
twist  disclination  solutions  [6,21]  for  linear-elastic  isotropic  bod¬ 
ies.  On  the  other  hand,  for  a  nonlinear  elastic  body  with  spatially 
constant  moduli,  components  of  the  integrated  symmetric  dis¬ 
placement  gradient  are  given  by  the  six  independent  integral  equa¬ 
tions 


CABCD  u(CJ})dV=-  CAECD  uBEuCJ}dV-  CBECD  uAEuCDdV 

J  y  J  y  J  y 

-  uE,cuEDdV 


1  qABCDEF 
’  2  ' 


I  uC,DuE,F^y 

J  v 


(39) 


Hence  from  Eqs.  (35)  and  (39),  components  of  F1  may  be  non- 
negligible: 


Fmn~  $mn~  V  1$mnab 


CAECD  f  uBEuCDdV 

J  v 


I  uA,EUC,1 

J  V 

L 


,dV+-CABCD 


+  iQBECD 

+  ^CABCDEF  |  uCDuEFdV 


j  uE,CuEJ 

J  V 


AV 


(40) 


Rank  4  elastic  compliance  SMNAB  in  Eq.  (40)  satisfies  [29,48] 


2CCZMfiS/lfiMA,  -  S CM(f>N+  (41) 

By  introducing  the  quantity 

qCDABEF  _  _  £ABCD  <pF  _  qABEF  fiOD  +  qAFCD  <pE  +  £BFCD g\E 
+  CEFADtfC  +  C  EFBD&C  +  C ABFD^E  +  QABEDgCF 

+  CABCDEF,  (42) 

relation  (39)  can  be  written  more  compactly  as  [44] 
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CABCDJ  Hc,D)dV=X- CABCDJ  UaEuEDdV-CABCDJ  uEEuCDdV 

_  I^cdabefJ  Uc  i)Ue  f(iv  (43) 

Correspondingly,  Eq.  (40)  can  be  replaced  with 

ABCDj  uCtEuEDdV 

-CABCDj  uEEuc,DdV-±C™ABEEJ  uCDuE  fdV 

(44) 

Notice  from  Eqs.  (40)  and  (44)  that  both  geometric  nonlinearity 
(quadratic  terms  in  displacement  gradients)  and  material  nonlin¬ 
earity  (third-order  elastic  constants)  contribute  to  F1,  and  that  F1 
does  not  necessarily  reduce  to  the  unit  tensor  when  the  third-order 
constants  vanish.  If  the  strain  energy  in  Eq.  (34)  is  extended  to 
incorporate  displacement  gradients  of  order  higher  than  3,  e.g., 
fourth-order  elastic  constants  [48],  then  the  effects  of  these  higher- 
order  terms  will  likewise  enter  the  right  sides  of  Eqs.  (40)  and 
(44).  When  elastic  moduli  are  not  spatially  constant,  for  example, 
in  a  multiphase  composite  or  body  with  foreign  inclusions,  Eqs. 
(40)  and  (44)  do  not  strictly  apply  since  in  that  case  elastic  coef¬ 
ficients  cannot  be  moved  outside  the  volume  integrals.  If  V  is  a 
polycrystal  with  randomly  oriented  grains,  then  uniform  isotropic 
elastic  properties  (i.e.,  effective  moduli)  can  be  assigned  to  each 
individual  crystal  as  an  approximation. 

Because  the  average  residual  elastic  deformation  results  from 
nonlinear  elastic  effects,  specifically  products  of  displacement 
gradients  of  order  2  in  Eqs.  (40)  and  (44),  the  contribution  of  the 
average  residual  elastic  deformation  F1  to  the  total  deformation 
gradient  F  of  Eq.  (15)  will  generally  be  small  in  conventional 
engineering  applications,  wherein  defect  densities  are  low  to  mod¬ 
erate.  For  example,  linear  elasticity  theory  is  generally  deemed 
valid  beyond  some  cutoff  distance  on  the  order  of  1-10  lattice 
parameters  from  the  dislocation  core  [29,57],  beyond  which  mag¬ 
nitudes  of  elastic  displacement  gradients  are  small  (less  than 
~0.1)  so  that  contributions  from  such  linear-elastic  regions  to 
terms  in  braces  in  Eq.  (44)  will  be  negligible.  On  the  other  hand, 
for  materials  with  very  large  dislocation  densities  in  which  core 
regions  comprise  a  substantial  fraction  of  the  volume  element, 
nonlinear  elastic  contributions  could  be  substantial,  in  which  case 
the  difference  F!-l  would  be  non-negligible.  The  contribution  of 
tensile  volumetric  and  deviatoric  parts  of  local  elastic  displace¬ 
ment  gradients  to  F1  would  also  be  limited  by  the  theoretical 
strength  of  the  crystal,  typically  on  the  order  of  10%  of  an  elastic 
shear  modulus  [57],  since  the  material  would  fracture  at  elastic 
strains  producing  tensile  or  shear  stresses  in  excess  of  the  theoret¬ 
ical  strength.  On  the  other  hand,  large  elastic  volumetric  compres¬ 
sions  are  usually  sustainable  in  crystals.  The  effect  of  residual 
elasticity  would  be  greatest  in  materials  whose  third-order  elastic 
constants  are  substantially  larger  than  second-order  constants;  in 
some  crystalline  solids,  representative  third-order  constants  can 
exceed  second-order  constants  by  an  order  of  magnitude  or  more 
[29], 

3.4  Average  Residual  Elastic  Volume  Changes.  Now  con¬ 
sider  the  average  volume  change  in  the  body  resulting  from  the 
field  of  local  residual  elastic  displacement  gradients  uA  B.  The  pre¬ 
ceding  derivations  (and  those  that  follow  in  Sec.  3.4)  apply  re¬ 
gardless  of  whether  or  not  the  average  shape  change  resulting 
from  uA  B  vanishes.  However,  for  a  crystal  containing  a  large  num¬ 
ber  of  randomly  oriented  defects  (e.g.,  dislocation  and  disclination 
lines  and  loops),  it  may  be  reasonable  to  assume  that  the  change  in 
dimensions  of  the  crystal  imparted  by  local  stress  fields  of  defects 


FmN  - 


exhibits  no  preferred  directions,  implying  that  the  crystal  under¬ 
goes  only  a  volume  change  and  no  shape  change.  The  average 
residual  elastic  volume  change  is  also  of  great  practical  interest 
because  it  can  be  easily  measured  experimentally  for  plastically 
deformed  crystals,  with  results  then  available  to  validate  the 
theory.  The  residual  elastic  deviatoric  shape  change,  on  the  other 
hand,  cannot  be  obtained  simply  from  the  deformed  shape  of  a 
sample  of  material  since  it  cannot  be  readily  delineated  from  the 
macroscopic  shape  change  resulting  from  dislocation  glide. 

The  change  in  a  differential  volume  of  a  body  in  coincident 
spatial  and  reference  coordinate  systems  is  measured  by  [44] 

dv/dV  =  det(<J^  +  uAB)  =  |[(x/|A)3  -  3xAAxBcxcB  +  2xBcxcAxAB\ 

=  1  +  u\  +  2  “  2  (45) 

where  terms  of  order  3  and  higher  in  displacement  gradients  are 
neglected  in  the  final  equality.  In  the  present  application  in  the 
context  of  Fig.  1,  du  is  the  volume  of  a  differential  subelement  of 
the  body  in  configuration  B,  and  dV  is  the  volume  of  that  subele¬ 
ment  in  B.  Integrating  Eq.  (45),  a  second-order  accurate  measure 
of  the  net  volume  change  attributed  to  residual  elastic  deformation 
within  the  crystal  is 

A V=J  uAAdV  +  |  [{uAA)2-uABuBA\dV  (46) 

The  right  sides  of  Eqs.  (40),  (44),  and  (46)  depend  on  local  re¬ 
sidual  elastic  displacement  gradient  fields  in  the  body.  For  crystals 
of  high  symmetry,  these  relations  can  be  further  reduced  by  ap¬ 
pealing  to  particular  forms  of  the  elastic  coefficients.  Specifically, 
for  cubic  crystal  systems  of  the  highest  symmetry  (Laue  group 
Cl),  second-order  moduli  consist  of  three  independent  coefficients 
and  the  third-order  moduli  of  six  independent  coefficients.  Con¬ 
sider  a  coordinate  system  with  axes  parallel  to  the  cube  axes  of  the 
crystal.  In  Voigt’s  notation  [58],  pairs  of  indices  11— >  1,  22— >2, 
33— >3,  23— >4,  13— >5,  12— >6,  and  2EAB^Ea(l  + SAB),  where 
Greek  indices  span  a=  1 ,2, ...  ,6,  nonzero  second-order  constants 
are 


C11-C22-C33,  C12-C13-C23,  C44  —  Cg5  —  C66  (47) 

and  nonzero  third-order  constants  are 

C111  =  C222  =  C333’  C144  =  C255  =  C366, 

Cll2=  Ctl3  =  C122  =  C223  =  C133  =  C233,  Q23,  (48) 

C155  =  C166  =  C244  =  C2 66  =  C344  =  C355’  C456 

The  following  notation  is  common  [11]: 

3B  —  C\i  +  2C\2i  2G=  Cn-  C12,  —2Gr  —  Cy\  — Ci2  —  2Caa 

(49) 

where  B  is  the  bulk  modulus  and  G  and  G'  are  the  shear  moduli. 
In  a  cubic  crystal,  the  first  term  on  the  right  of  Eq.  (34)  can  be 
written  as  [44] 

W=\Cabc°ua,buc„=\B{uaaY 

+  G[(z/>fi4  +  «>/)/ 2  -  (z/A)2/ 3] 

+  G'[(uabuba  +  uabuab)/2  -  (ulA)2  -  (1 r2)2  -  («j)2] 

=  WD+WS+W'S  (50) 

where  Wo  is  the  strain  energy  of  dilatation,  and  VF;  and  W's  result 
from  shape  changes  of  the  crystal.  In  an  isotropic  material  (e.g.,  a 
polycrystal  with  no  preferred  orientations),  the  number  of  con¬ 
stants  in  Eq.  (47)  is  further  reduced  to  2  according  to 

2C44  =  C11-C12  (51) 

and  in  Eq.  (48)  reduced  to  3  according  to 
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2C\a4  —  C 


144-  ^  112  —  C123 


4C,5s  =  C 


155  —  '-111-  Cil2 


8C456-  Cm  -  3C112  +  2C123 


(52) 


FdV 


(53) 


WdV 


Ca™  u(CD)dV=  -3  1 

J  v  J  V 

'  fA.CD  I  /  E  £  341/ 

J  \UC,EU.,D~  uC,Du.,E)dV 

1  AA.CDEF  I 

2'  J „ 


+  -l 
2 


uc,due,fl iV 


(54) 


For  a  cubic  crystal  satisfying  Eqs.  (47)-(50),  this  reduces  to 

3 bJ  uDDdV=-  3J  WdV  -  J  [(uAA)2  -  uAguBA]dV 

J  uC,DuE,FdV 

J  v 


_AA.CDEF 
'2  A 


(55) 


AV=-~  WdV-—CAl 
Bjv  6 B  A 


pf 

I  UC,DUE„ 

J  w 


FdV  (56) 


Letting  det(x°A)  =  u,  where  u>0,  the  rightmost  term  of  Eq.  (56) 
can  be  written  as  [44] 


1 

AV=- 

B 


SB 

—  -  1 

dp 


l 


B  dG 


WDdV  +  I  —  —  -  1 
v  VG  dp 


L 


WsdV 


,  B  dG' 

+  — - -1 

kG'  dp 


W'sdV 

J  V 


(58) 


WD  =  V~ll  WDdV ,  ws=\r'l  WsdV,  w's=\r'l  W’sdV 
J  v  J  v  J  V 


Thus  in  an  isotropic  material,  B  and  G  of  Eq.  (49),  and  WD  and  Ws 
of  Eq.  (50)  are  unchanged,  G'=0  in  Eq.  (49),  and  VT^=0  in  Eq. 
(50).  Summing  over  the  first  two  indices  of  Eq.  (39),  and  using  the 
first  of  Eq.  (50)  gives 

C aCDJ  u(CJ))dV=-4^  WdV- fa™  f  uECuEDdV 

'  ,'A.CDEF  I 

~  7  J  UC,DUEJ 

Appealing  to  Eq.  (43), 


(59) 


and  the  quantity 


J=\+\VIV  (60) 

the  normalized  volume  change  for  a  self-equilibrated,  cubic  non¬ 
linear  elastic  solid  is  found  as 


1 

:1  +  - 


dB  \-  l  B  dG  \-  I  B  dG’ 

~  -1)^0+  -1  Ws+(— — -1|W( 


B\_\dp  ! 


G  dp 


a  dp 


(61) 

If  the  cubic  second-order  elastic  coefficients,  their  pressure  de¬ 
rivatives,  and  each  of  the  average  strain  energy  densities  of  Eq. 
(59)  are  known,  normalized  residual  volumetric  deformation  can 
be  computed  from  Eq.  (61).  If  strain  energy  densities  are  all  posi¬ 
tive,  and  if  coefficients  of  the  energy  densities  in  Eq.  (61)  are  all 
of  the  same  sign  (e.g.,  positive),  then  the  overall  volume  change 
will  be  of  that  sign  (e.g.,  positive).  If  only  the  total  strain  energy 
density  W=WD+WS+W'S  is  known,  for  example,  from  experi¬ 
ments  [11,45],  then  Eqs.  (50)  and  (61)  can  be  combined  to  estab¬ 
lish  bounds  on  the  normalized  volume  change  [44].  Such  bounds 
have  been  validated  for  several  polycrystalline  cubic  metals  [11], 
wherein  the  volume  changes  were  always  found  positive  by 
theory  and  experiment.  For  an  isotropic  solid,  Eq.  (61)  reduces  to 


J=l  +  -[  —  ~l\WD+-\—--\Ws 
B  \  dp  )  D  G\ dp  Bj  s 


j JdG  G 


(62) 


In  an  isotropic  body,  pressure  derivatives  of  tangent  bulk  and 
shear  moduli  B  and  G  in  the  undistorted  reference  configuration 
are  related  to  third-order  elastic  constants  by  [29] 


Then  from  Eq.  (46),  the  net  volume  change  for  the  element  of 
reference  volume  V  is 


dB 
-B  — 

dp 

dG 

-B  — 

dp 


P= 0 


8 

:  C123  +  2C144  +  -C456 


p= 0 


1  4 

:  B  +  —G  +  C144  +  — C456 


(63) 


I  UC,DUE,F dV 

J  V 

v  dB  f  v  dG  f  V  dG'  f  , 

= - WDdV+ - WsdV+— -  W'sdV 

B  dvj  y  GdvJv  G'  dvjy  5 

(57) 

Setting  dl dv=-{BI  v)  d / dp,  where  p  is  the  Cauchy  pressure,  and 
substituting  Eqs.  (50)  and  (57)  into  Eq.  (56), 


where  moduli  B,  G,  and  G'  and  their  derivatives  with  respect  to 
pressure  are  all  evaluated  at  a  stress-free  reference  state.  After 
defining  average  strain  energies  on  a  per-reference-volume  basis 
as 


Formula  (61)  for  cubic  crystals  is  attributed  to  Toupin  and  Riv- 
lin  [44],  while  Eq.  (62)  was  developed  earlier  by  Zener  [8],  Iso¬ 
tropic  formula  (62)  for  the  case  when  the  dilatational  energy  van¬ 
ishes  agrees  with  that  of  Holder  and  Granato  [10]  obtained  by 
simple  thermodynamic  arguments: 

dg  1  (dGe  Ge\  1  I  dGe  Ge\  ,  N 

AVIV=  —  =  —  — g  =  —  — --  —  )EPt  (64) 
dp  Ge\  dp  B  r  Ge\  dp  bJh1 

where  g  is  the  Gibbs  free  energy  change  from  defects,  per  unit 
reference  volume  that  depends  on  pressure  p  and  temperature  9 , 
and  Ge  is  an  effective  elastic  constant  that  depends  on  the  math¬ 
ematical  form  of  the  strain  energy  of  the  particular  defect.  In  the 
final  term  of  Eq.  (64),  E  is  the  elastic  strain  energy  per  unit  length 
of  the  defect,  and  pT  is  the  total  length  per  unit  volume  of  the 
defect.  For  a  number  of  metallic  crystals,  Holder  and  Granato  [10] 
found  that  delineation  of  dilatational  and  deviatoric  energies  and 
effects  of  anisotropy  had  little  effect  on  volume  changes  from 
straight  dislocations  predicted  using  Eq.  (64). 

The  source  of  the  local  displacement  gradient  and  residual 
strain  energy  in  the  body  to  this  point  has  been  arbitrary,  so  long 
as  Eqs.  (20),  (23),  and  (34)  apply;  that  is,  traction  is  continuous 
across  all  internal  surfaces,  the  body  is  self-equilibrated  such  that 
external  traction  is  absent,  and  the  crystal’s  local  constitutive  re¬ 
sponse  is  described  by  hyperelasticity  with  terms  of  order  higher 
than  3  in  the  displacement  gradients  neglected  in  the  strain  energy. 
Defects  induce  such  displacement  gradients  and  residual  strain 
energies.  These  may  include,  for  example,  dislocation  and  discli- 
nation  lines  and  loops,  stacking  faults,  grain  boundaries,  twin 
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boundaries,  and  slip  bands.  On  the  other  hand,  volume  changes 
attributed  to  point  defects  (e.g.,  substitutional  atoms,  interstitials, 
and  vacancies),  phase  transformations,  voids,  and  open  cracks  are 
not  considered  in  the  present  treatment,  since  dimensional 
changes  computed  according  to  the  present  theory  account  only 
for  volume  and  shape  changes  resulting  from  stress  fields  of  de¬ 
fects  and  not  volume  and  shape  changes  associated  with  defects 
themselves.  However,  the  additional  volume  change  induced  by  a 
point  defect — the  volume  change  in  addition  to  the  misfit  dilata¬ 
tion  in  a  sphere-in-hole  model — attributed  to  elastic  nonlinearity 
can  be  estimated  from  Eq.  (62)  for  an  isotropic  elastic  body  [19]. 
Anisotropy  cannot  be  directly  addressed  in  Eq.  (61)  to  describe 
the  effects  of  defects  wherein  moduli  may  change  with  position, 
for  example,  grain  boundaries  and  twin  boundaries  across  which 
lattice  orientations  may  differ.  However,  isotropic  approximation 
(62)  could  be  used  as  an  estimate  in  these  cases,  as  in  previous 
applications  toward  polycrystals  [43]  and  single  crystals  of  lower 
symmetry  [15]. 

Also  neglected  in  the  foregoing  continuum  elastic  analysis  are 
explicit  effects  of  defect  cores  on  residual  deformation.  For  ex¬ 
ample,  elastic  strain,  stress,  and  strain  energies  imparted  by  Volt- 
erra  line  defects  diverge  as  the  radial  distance  from  the  line 
shrinks  to  zero  [6],  and  even  nonlinear  elasticity  is  inadequate  for 
describing  the  nonconvex  energy  distribution  imparted  by  the 
highly  distorted  crystal  structure  within  defect  cores.  In  the 
present  treatment,  one  can  imagine  each  defect  line  to  be  sur¬ 
rounded  by  a  traction-free  cylindrical  boundary  delineating  the 
core  volume  from  the  surrounding  elastic  continuum  of  reference 
volume  T.  In  linear  elasticity,  stresses  and  strains  resulting  from 
traction  acting  on  the  cylindrical  boundary  decay  as  R~ 3,  where  R 
is  the  distance  from  the  defect  line,  and  hence  are  usually  ne¬ 
glected  [29],  The  elastic  continuum  is  self-equilibrated  in  this  ap¬ 
proximation,  and  Eq.  (36)  applies. 


AT 

L 


1  /  SB 


edge 


=  -—-l£*  + 


B\dp 


]_(dG  G 
G  \  dp  ~  B 


l-v-2  iP-ldB  \  2-2v+2v2/dG 


3  L  (1  -v)B  \dp 


(1  -  v)G  \  dp 


where  Ee=E^+Ee  is  the  total  elastic  line  energy  of  an  edge 
dislocation. 

Analogously,  as  denoted  by  ES=E-E,  the  strain  energy  per  unit 
line  length  L  of  straight,  noninteracting  screw  dislocations,  parti¬ 
tioned  as  [43] 


Letting 


4=o, 


(69) 


WSV=ESSL  =  ESL  (70) 

where  Es  is  the  total  elastic  line  energy  of  a  screw  dislocation,  and 
multiplying  Eq.  (62)  by  VI L,  the  volume  change  per  unit  length  of 
straight  screw  dislocation  lines  is  [29,43] 


AT  G\ 

L  screw  G  V  dp  Bl 


(71) 


Now  let  the  total  dislocation  line  density  per  unit  reference 
volume  be  given  by  the  sum 


Pt  =  L/V  =  pE  +  ps  (72) 

where  now  L  is  the  total  length  of  edge  and  screw  dislocations  in 
volume  T,  pE  is  the  density  of  edge  dislocations,  and  ps  is  the 
density  of  screw  dislocations.  Define  from  Eq.  (65) 


3.5  Straight  Edge  and  Screw  Dislocations.  Consider  volume 
changes  imparted  by  dislocation  lines  in  the  isotropic  approxima¬ 
tion.  The  energy  per  unit  length  (linear-elastic  energy  plus  core 
and  interaction  energies)  of  a  straight  dislocation  in  an  infinite 
isotropic  medium  is  written 


1  -  v  (edge  dislocation) 

1  (screw  dislocation) 

(65) 

where  b  is  the  magnitude  of  the  Burgers  vector,  v=(3B 
-2G)I(6B  +  2G)  is  Poisson’s  ratio,  R  is  the  radial  distance  from 

the  dislocation  core,  Rc  is  the  radius  of  the  dislocation  core,  and  E 
is  a  correction  that  accounts  for  core  energy,  line  curvature,  inter¬ 
action  energies  from  other  defects  and  boundaries,  and  stacking 
faults  associated  with  partial  dislocations. 

Denoted  by  Ee=E-E  is  the  elastic  strain  energy  density  per 
unit  line  length  L  of  straight,  noninteracting  edge  dislocations, 
separated  into  dilatational  and  shear  contributions  as  [43] 


E  =  - 


Gb2 


4ttK  \R, 


ln|* 


1 c 


Gb2  \-v-2v1  I  R\ 
127 T  ( l-v f  \RC)’ 


F  Gb-2-2v+2ir  /  R 

Es  = - : - ry — In  — 

1277  ( l-v )2  \RC 

(66) 


Letting 


WdV  =  EedL,  WsV=EesL  (67) 

and  multiplying  Eq.  (62)  by  VI L,  the  volume  change  per  unit 
length  of  edge  dislocation  line  is  [29,43] 


Et=KEe  =  Es  = 


(73) 


and  let  x=Pe!  Pt  be  the  fraction  of  pure  edge  dislocations  in  the 
total  dislocation  density.  Supeiposing  volume  changes  from  Eqs. 
(68)  and  (71),  and  using  Eq.  (72),  Eq.  (62)  becomes 
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l  +  i^ 


+  Ps 

edge 
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SA'L  (1  -v)B  \dp 
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2-2v+2  lrldG  G 
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(74) 


Relation  (74)  yields  the  volume  change  associated  with  straight, 
noninteracting  pure  screw  and  edge  dislocation  lines  in  an  isotro¬ 
pic  body  with  homogeneous  elastic  properties.  While  terms  in 
braces  in  Eq.  (74)  include  nonlinear  elastic  effects  (i.e.,  pressure 
derivatives  of  the  elastic  moduli  or  third-order  elastic  constants 
from  Eq.  (63)),  dislocation  line  energy  E T  does  not  account  ex¬ 
plicitly  for  nonlinearity  or  energy  of  dislocation  cores.  Thus,  Eq. 
(74)  represents  the  product  of  third-order  nonlinear  elasticity 
(terms  in  braces)  and  linear  elasticity  (defect  energy).  Products  of 
purely  linear-elastic  origin  would  yield  null  volume  change,  as 
implied  by  Eq.  (38),  while  higher-order  products  of  nonlinear  ori¬ 
gin  (e.g.,  products  of  pressure  derivatives  of  moduli  with  contri¬ 
butions  of  nonlinear  elasticity  to  dislocation  energies)  are  ne¬ 
glected.  Relation  (74)  can  be  used  directly  to  compute  the  volume 
change  in  a  (poly)crystal  if  line  densities  of  edge  and  screw  dis¬ 
locations,  their  energies  per  unit  length,  and  the  requisite  elastic 
constants  are  known.  The  preceding  analysis  assumes  that  defect 
densities  are  negligible  in  the  reference  configuration.  When  dis¬ 
location  densities  are  substantial  in  the  initial  state,  the  predicted 
volume  change  between  the  initial  and  final  configurations  corre¬ 
sponds  only  to  the  change  in  defect  density  that  occurs  during 
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deformation,  rather  than  that  associated  with  the  absolute  density 
of  dislocations. 

Dislocation  line  densities  in  the  preceding  developments  are 
defined  per  unit  reference  volume  (equivalent  in  configurations  B0 
and  B  of  Fig.  1).  When  the  dislocation  line  density  pT  is  measured 
in  the  externally  unloaded  but  internally  stressed  configuration  B, 


then  pT^>JpT  in  Eq.  (74)  such  that 


1  -  v-2v r!  dB 
3Kl  (1  -v)B  l  dp 
1 ~x(dG  G) 

G  [dp 


dB  \  2-2v  +  2i?ldG 


(1  -  v)G  \  dp 


ETpT 


However,  the  distinction  between  configurations  used  to  define  the 
dislocation  line  density  will  have  a  trivial  effect  on  computed 
values  of  J  for  small  volume  changes  (1+AWV~1/(1-AW  VO), 
and  in  these  cases  preferential  use  of  Eq.  (75)  over  Eq.  (74)  is 
probably  not  warranted  given  the  uncertainty  in  dislocation  den¬ 
sities  that  can  be  measured  experimentally.  The  dislocation  line 
energy  in  Eq.  (73)  depends  on  the  core  radius  Rc  and  the  cutoff 
radius  R.  A  typical  approximation  is  [57] 


ET=aGb2 


(76) 


with  0.5:£a'=ln(F/Rc)/4'7r£  1.0.  Introducing  the  dimensionless 
quantities 


ETpT 


1  -  v-2ir(  dB 


3K  L  (1  -  v)B  \dp 


dB  \  2-2v+2v1/dG 


(1  -  v)G  \  dp 


_EIp1(dG  G 
G  [  dp  ~  B 

the  volume  change  in  Eq.  (74)  becomes 


\VIV=J-1=AeX  +  As(1-x)  (78) 

Thus,  A  VI V  for  a  mixture  of  straight,  noninteracting  edge  and 
screw  dislocations  will  fall  between  the  volume  change  AE  result¬ 
ing  from  the  same  density  of  pure  edge  dislocations  and  the  vol¬ 
ume  change  As  resulting  from  the  same  density  of  pure  screw 
dislocations. 


3.6  Approximate  Volume  Changes.  Returning  to  the  overall 
dimensional  changes  of  the  body  induced  by  averaged  local  dis¬ 
placement  gradients,  consider  Eq.  (35).  Defining  the  quantities 

K=\Tl\  uAAdV,  Z2=W2  uABdV  I  uBAdV  (79) 
J  v  J  v  J  V 

and  neglecting  terms  of  orders  3  and  higher  in  displacement  gra¬ 
dients,  the  analog  of  Eq.  (45)  is 

det(FAS)  =  1  +K  +  K2/ 2  -  L2/ 2  (80) 

Recall  from  Eqs.  (46)  and  (60)  that 

7=  1  +  K+  i  J  [(i/A)2  -  uABuBA]dV  (81) 

Thus, 

7  =  det(F^)  -  1-{K2  -  L2)  +  i  J  l(uAA)2  -  uABuBA\dV  (82) 

so  that  at  least  to  first-order  in  displacement  gradients,  J 
=  det (Fab).  Possible  differences  between  7  and  det(FAB)  may  arise 
because  the  determinant  operation  and  volume  integration  do  not 
commute.  When  products  of  displacement  gradients  and  their  av¬ 
erages  are  small,  F7A-<?^  =  F~det(FAB)-l  and  7~det(FAB).  Re¬ 


call  that  EAB=FAB)  is  a  stretch  with  no  rotation.  Even  though  in 
some  cases  defects  impart  no  local  dilatation  in  the  linear-elastic 
approximation  (e.g.,  Volterra  screw  dislocations  in  isotropic  sol¬ 
ids),  a  volume  change  can  still  result  (e.g.,  Eq.  (71))  because  of 
nonlinear  elastic  effects. 

3.7  Examples.  Table  1  presents  strain  energies  per  unit  line 
length  E  of  a  number  of  defects:  straight  edge  and  screw  disloca¬ 
tions  discussed  already  in  Sec.  3.5;  circular  screw  and  prismatic 
dislocation  loops  of  radius  R  and  core  radius  Rc,  and  straight 
wedge  disclinations  and  circular  twist  and  wedge  disclination 
loops  of  Frank  vector  to,  radius  r,  and  core  radius  rc.  Energies  all 
correspond  to  defects  in  infinite,  isotropic  linear-elastic  solids. 
Analytical  formulas  for  energies  follow  from  references  quoted  in 
Table  1,  while  normalized  volume  changes  (AWV)/(pjd2)  com¬ 
puted  for  copper  follow  from  Eq.  (68)  for  edge  dislocations  and 
from  Eq.  (62)  for  the  other  defects,  using  material  properties  listed 
in  Table  2.  For  defects  besides  edge  dislocations,  strain  energy  is 
treated  as  deviatoric,  leading  to  the  linear  relationship  between 
volume  change  and  defect  density  given  in  Eq.  (64).  This  is  a 
rigorous  assertion  for  screw  dislocation  loops  and  twist  disclina¬ 
tion  loops  but  not  prismatic  loops  or  wedge  disclinations,  though 
it  has  been  used  elsewhere  for  these  defects  [18,61]  and  should 
provide  a  reasonable  approximation  to  the  volume  change  for 
many  metals  [10],  including  Cu.  Again,  pT  is  the  line  length  per 
unit  reference  volume  of  the  defect,  and  a  is  the  lattice  spacing  of 
the  conventional  unit  cell  in  a  perfect  crystal  at  room  temperature. 
Experimental  results  for  polycrystalline  Cu  are  provided  for  com¬ 
parison  [45];  the  range  of  volume  changes  reported  in  Table  1  for 
the  experiments  corresponds  to  estimates  of  dislocation  line  length 
from  energy  measurements  after  plastic  compressive  strains  rang¬ 
ing  from  0.3  to  0.7.  The  nature  of  the  defects  (e.g.,  edge  versus 
screw  and  straight  lines  versus  loops)  was  not  reported  in  that 
experimental  investigation. 

From  Table  1,  the  normalized  volume  changes  predicted  for 
various  kinds  of  dislocations  agree  with  the  experimental  findings 
within  a  factor  of  ~2  and  are  very  close  for  edge  dislocations  and 
prismatic  loops.  Volume  changes  per  unit  defect  length  in  Cu  are 
positive  (i.e.,  defects  cause  expansion)  and  small.  For  example,  a 
1%  volume  increase  would  require  an  immense  density  of  edge 
dislocations  on  the  order  of  pT~0.0la~2~  1017  m~2,  correspond¬ 
ing  to  an  average  dislocation  spacing  on  the  order  of  10a.  For 
small  volume  changes  (AWV<1%),  det  F1  should  be  expected 

provide  an  accurate  approximation  of  J  according  to  Eq.  (82). 
Disclination  energies  and  volume  changes  are  comparable  to  those 
for  dislocations  for  the  small  cutoff  radii  or  small  loops  consid¬ 
ered  here  (?■=  10rc)  but  would  diverge  quickly  at  large  distances  r 
from  the  core  of  straight  disclinations  or  for  larger  disclination 
loops.  Elastic  anisotropies  of  individual  grains,  defect  core  ener¬ 
gies,  and  interaction  energies  among  defects  and  between  defects 
and  the  external  boundary  of  the  body  (e.g.,  image  forces)  are 
neglected  in  this  application  of  the  model;  such  effects  presum¬ 
ably  contribute  to  discrepancies  between  theory  and  experiment. 

The  results  of  Table  1,  along  with  experimental  data  summa¬ 
rized  by  Wright  [11],  suggest  that  volume  changes  resulting  from 
residual  elasticity  associated  with  dislocations  should  be  small  in 
metallic  crystals  of  cubic  symmetry  deformed  in  compression  or 
shear  to  strains  on  the  order  of  unity,  under  quasistatic  conditions. 
The  experimental  results  in  Table  1  [45]  correspond  to  dislocation 
densities  on  the  order  of  pT~  1015  m-2.  leading  to  volume 
changes  on  the  order  of  AW  V~  10-4.  Such  small  volume  changes 
would  seem  inconsequential  in  the  context  of  measured  yield 
properties  of  materials  in  unconfined  loading,  when  specimens  are 
free  to  expand  or  contract  laterally.  However,  under  lateral  con¬ 
finement,  for  example,  uniaxial  strain  conditions  occurring  in 
shock  loading  [14,15],  small  volume  changes  can  drastically  af¬ 
fect  the  measured  hydrostatic  pressure.  Volume  changes  cannot  be 
accommodated  by  dislocation  glide,  as  noted  in  Sec.  2.1;  hence, 
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Table  1  Defect  energies  and  residual  volume  changes  in  copper 


Defect 

Energy/length  E 
(nJ/m) 

Volume  change 
(AV 1 V)  /  (pja2) 

Edge  dislocation3: 

Gb 2  R 

E-  In 

4tt(1  -  v)  Rc 

5.29 

0.89 

Screw  dislocation3: 

Gb2  R 

E  = - In — 

4t r  Rc 

3.39 

0.51 

Screw  dislocati 
Gb 2 

E~ - 

127t(1  -  v ) 

an  loopb: 

8  R 

3(1  -  i/)Ln —  -  (9  -  lOn) 
Rc 

3.21 

0.49 

Prismatic  dislocation  loopb: 

5.32 

0.81 

Gb2 


47t(1  -  v) 


8  R 

In— -2 
Rr 


Wedge  disclination': 


E=- 


Gu2r 


16  tt(  1  -  v) 


1- 


-  -4  -  1-  - 


In2  — 


4.24 


0.64 


Twist  disclination  loop1": 


Gurr 2 


47 r 


8r  8 
In—  -  - 

rr  3 


15.8 


2.4 


Wedge  disclination  loop' 

Garr 
E 


8  r  8 

In - 

3 


12.4 


1.9 


Experiment' 


0.8-1. 1 


'Reference  [21]. 
bReference  [59]. 
'Reference  [60]. 
dReferences  [61,62]. 
'Reference  [45]. 


any  residual  inelastic  volume  increase  must  be  compensated  by  an 
elastic  volume  decrease  and  corresponding  hydrostatic  pressure. 
Extremely  high  dislocation  densities,  pT~  1016  m-2,  have  been 
observed  for  metals  such  as  Cu  deformed  in  shock  conditions 
[14],  A  dislocation  density  of  1016  m-2  would  lead  to  an  expan¬ 
sion  on  the  order  of  0.1%,  which  would  in  turn  require  an  offset¬ 
ting  pressure  on  the  order  of  0.1%  of  the  bulk  modulus,  or 
— 0. 1 5  GPa  in  Cu.  Such  contributions  of  defects  to  pressure  or 
volume  change  would  manifest  implicitly  in  the  measured  equa¬ 
tion  of  state  (i.e.,  pressure-volume-temperature  relationship),  bulk 
modulus,  and  pressure  derivatives  of  the  bulk  modulus  obtained 
from  high-pressure  experiments.  Substantial  effects  of  large  dislo¬ 
cation  densities  on  tangent  elastic  moduli  of  metallic  single  crys¬ 
tals  deformed  in  uniaxial  strain  have  been  predicted  using  atom¬ 
istic  methods  [51],  Very  large  dislocation  densities  comprising 
dislocation  walls  in  subdivided  grains  have  been  observed  in  met¬ 
als  subject  to  severe  plastic  deformation  [63];  residual  elastic  vol¬ 
ume  changes  could  be  significant  in  highly  defective  regions  of 
such  materials. 

While  closed-form  analytical  solutions  are  not  available  for  the 


deviatoric  part  of  F1,  magnitudes  of  its  deviatoric  components  are 
expected  to  be  on  the  same  order  as  its  volumetric  part,  i.e.,  small 
except  in  cases  where  defect  densities  are  extremely  large.  Meant 
by  small  are  conditions  Fab- Sab  <0.001  for  defect  densities  pT 
<  1015  m-2  in  cubic  metals  such  as  copper,  though  larger  residual 
elastic  deformations  would  be  expected  in  crystals  that  feature 
stronger  pressure  sensitivity  of  the  elastic  coefficients.  Because 
deviatoric  deformation  can  be  accommodated  by  plastic  slip,  there 
seems  to  be  no  capacity  available  to  delineate  contributions  from 
residual  elasticity  versus  dislocation  flux  to  the  plastic  properties 
of  materials  (e.g.,  yield  and  flow  properties)  measured  in  experi¬ 
ments  on  single  or  polycrystals.  The  effect  of  F1  is  implicitly 
included  in  such  experimental  measurements.  The  effects  of  F1 
could  in  principle  be  delineated  via  numerical  simulations  incor¬ 
porating  atomistic  or  discrete  dislocation  dynamics.  Presumably, 
contributions  of  residual  elasticity  to  overall  deformation  gradient 
F  are  small  in  conventional  engineering  applications  wherein  de¬ 
fect  densities  are  small  to  moderate,  since  omission  of  F1  in  stan¬ 
dard  crystal  plasticity  models  [53]  does  not  seem  to  inhibit  accu- 


Table  2  Properties  for  copper3 


G 

(GPa) 

B 

(GPa)  v 

dG/dp 

dB/dp 

a 

(nm) 

b 

(nm) 

r 

(nm) 

Rc 

co  (nm) 

R/ Rc 

rtrc 

47 

152  0.36 

0.8 

4.4 

0.362 

0.256 

1.0 

tt/2  1.0 

106 

10 

dReference  [43]  (elastic  coefficients  and  their  pressure  derivatives). 
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racy  of  such  models.  On  the  other  hand,  interacting  stress  fields  of 
dislocations  strongly  affect  yield  and  strain  hardening  [57];  fur¬ 
thermore,  spatial  gradients  of  F1  need  not  always  be  small  even  if 
F1  itself  is  small  relative  to  the  unit  tensor.  Gradients  of  residual 
elastic  deformation,  in  addition  to  gradients  of  lattice  rotation 
[30],  can  be  associated  with  the  density  of  geometrically  neces¬ 
sary  dislocations  and  kinematic  hardening  in  metals  [34,35]. 


tions),  the  average  virial  stress  for  a  group  of  atoms  is  often  con¬ 
sidered  equivalent  to  the  average  Cauchy  stress  [46]: 


«•=  2^2  r(u)  <8  f(l;>  (90) 


4  Atomistic  Approach 

The  volume  element  of  the  crystalline  solid  is  now  treated  as  a 
collection  C  of  discrete  atoms  rather  than  a  continuum.  Prior  to 
discussion  in  Sec.  4.4  of  the  self-equilibrated,  relaxed  intermedi¬ 
ate  configuration  and  means  to  compute  F1  entering  Eq.  (17) 
atomically,  relevant  background  discussion  pertinent  to  lattice 
statics,  atomic  stress  measures,  and  atomic  interactions  are  pro¬ 
vided  in  Secs.  4. 1-4.3. 


4.1  Lattice  Statics.  The  atoms  occupy  a  perfect  lattice  in  the 
volume  element  corresponding  to  configuration  B0  of  Fig.  1 .  As 
noted  already  in  Eq.  (1),  the  position  of  each  atom  i  in  configu¬ 
ration  B0  is  expressed  by  vector  R^,  with  angular  brackets  used  to 
denote  the  atomic  labels.  The  position  vector  of  atom  i  in  current 
configuration  B  is  r^\.  In  classical  Newtonian  mechanics  de¬ 
scribed  via  rectangular  Cartesian  coordinates,  the  balance  of  linear 
momentum  for  particle  i  is 

(83) 

where  is  the  mass  of  atom  i,  are  the  external  forces  that 
may  be  explicitly  time  dependent,  and  are  the  conservative 
forces  satisfying 


b^_ 
n  b 
dr«) 


(84) 


Potential  energy  <I>  is  a  function  of  all  atomic  coordinates: 


^  =  3>(r<i))ie£  (85) 

In  the  context  of  lattice  statics,  external  and  internal  forces  bal¬ 
ance.  In  the  absence  of  external  forces, 


dr\ 


=  0 


<•> 


(86) 


meaning  that  the  vector  sum  of  all  interatomic  forces  acting  on 
atom  i  is  zero.  A  more  specific  form  of  the  potential  that  maintains 
translational  invariance  of  the  energy  depends  only  on  the  relative 
separation  between  each  pair  of  atoms  r(y)=r(/)_r(;): 

<T>  =  tlJ(r<i>)‘je'C  (87) 

i*i 

leading  to 

<94> 

4>=*V-  (88) 

dr«j) 

with  components  of  an  interaction  force  between  atoms  i  and 
j.  In  general,  Eq.  (87)  accounts  for  three-body  and  higher-order 
interactions,  for  example,  angles  subtended  by  vectors  r^)  among 
interacting  sets  of  three  or  more  atoms.  For  materials  featuring 
central  force  interactions  only,  Eqs.  (87)  and  (88)  reduce  to 


i,j  e  C 

i*j 


r(ij)  =  (rW>  ’  r<i/>)1/2’ 


%>  - 


d<\>  r{ij) 
dr<a)  rdj ) 


(89) 


4.2  Atomic  Stress  Measures.  A  number  definitions  of  stress 
based  on  atomic  quantities  have  been  suggested  [46,50,64],  In  the 
context  of  lattice  statics  (i.e.,  null  atomic  velocities  and  accelera- 


where  1 1  is  the  total  volume  occupied  by  the  group  C  of  atoms, 
i.e.,  the  sum  of  all  atomic  volumes,  in  spatial  configuration  B.  The 
summation  in  Eq.  (90)  is,  in  fact,  a  double  sum,  proceeding  over 
all  pairs  of  distinct  (i  A  j)  atoms.  The  factor  of  2  arises  in  Eq.  (90) 
because  in  the  notation  used  here,  summation  proceeds  over  all 
i  A  j,  *(ij)=-T(ji),  and  f{ij)=~f(jiy  F°r  systems  with  central  force 
interactions  only,  symmetry  of  the  stress  tensor  follows  trivially 
from  Eqs.  (89)  and  (90): 


^2- 


(94> 


m  i*j  rdj > 


dr. 


\ij ) 


r(ijY(U)  - 


20 


-2 


1 

r<ij)dr(ij) 


rb  ra  _  J» 


(91) 


Now  consider  a  homogeneous  deformation  F  applied  to  all  atoms 
in  C ,  leading  to  r^  =  FR^  [22,38],  with  R^=R^-Rq-^  fixed 
separations  between  atoms  of  the  reference  lattice.  Rigid  body 
translations  are  omitted  in  this  homogeneous  deformation.  Letting 
the  average  strain  energy  per  unit  reference  volume  of  the  lattice 
'F  satisfy  flo^(F)  =  <&(r^-)(F,R^-))),  the  analog  of  Eq.  (33)  is 


paA  _  _  J_y  //<l*  dr<‘f>  _  _J_  y  gih  j.A 

dFaA  2O0fj  d^ij}dFaA  20^  drb{ij)  {ij) 

(92) 

where  O0  is  the  total  volume  occupied  by  the  atomic  system  in  B0. 
(The  atomic  volume  in  a  perfect  lattice  is  equal  to  Oq  divided  by 
the  number  of  atoms  comprising  £.)  Relation  (92)  is  a  conse¬ 
quence  of  the  affine  deformation  of  all  atoms  that  results  in 
dr^ij)/ dF aA=  SP-bR^ij).  The  average  atomic  Cauchy  stress  or  average 
static  virial  stress  is  then 


=  det(yi1A)RMF()i 


a* 


<?<i> 

drc,:! 


m 


1  ^ 

III  dr(ij) 


(93) 


in  agreement  with  Eq.  (90).  The  present  analysis  applies  only  to 
materials  described  by  potentials  of  the  general  form  (87)  and  for 
which  stress  relations  (90),  (92),  and  (93)  are  appropriate.  For 
example,  metals  that  can  be  modeled  by  combinations  of  pair 
potentials  and  multibody  potentials  such  as  the  embedded  atom 
method  are  included.  Not  admitted  are  piezoelectric  crystals  (e.g., 
noncentrosymmetric  ionic  solids),  some  of  whose  atoms  (sublat¬ 
tices)  may  display  a  relative  shift  when  polarized,  and  for  which 
hyperelastic  relation  (92)  may  not  be  sufficient.  In  these  cases, 
atomic  vibrations  may  provide  insight  into  origins  of  stress  and 
material  coefficients  [22,38,65].  Also  excluded  from  the  analysis 
are  noncentrosymmetric  polyatomic  lattices  such  as  diamond  and 
silicon  that  may  incur  inner  displacements  [54], 

4.3  Harmonic  and  Anharmonic  Interactions.  Expanding 
the  potential  energy  of  Eq.  (87)  in  a  series  about  the  reference 
state,  let 
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,  ,  '7‘f1 

®  =  <i,o+2j  tjt 


s*j 


%> 


» A  + 


is 


a1® 


i*j  dr0WdrU 


a  b 

%})Ci(kl) 


f{ij) = E  4*44  ^ = ^E  ^KkiA^U 


k+l 


i*j 

k+l 


4e 


31  "  dra{ij)drb{kl)dr\mn) 


9tiJ)9(kt)9{mn)  +  •  ■  ■ 


(94) 


*#/ 

m^n 


paA=^KkM^j) 

Zll0 i*j 


(100) 


where  <t>0  is  the  cohesive  or  ground  state  energy  of  the  reference 
lattice,  q(;i) =!•(;;•> -R(;y)  is  the  relative  displacement  between 
atomic  pairs,  and  the  second  term  vanishes  by  Eqs.  (86)  and  (88) 
since  the  reference  configuration  B0  is  explicitly  chosen  free  of 
external  and  internal  forces: 


<?<!> 

^4> 


=  0 


(95) 


Following  the  scheme  first  used  in  Eq.  (90),  repeated  use  of  the 
summation  symbol  is  omitted  in  Eq.  (94).  Hence,  summation  in 
the  second  term  on  the  right  side  of  Eq.  (94)  applies  over  two  sets 
of  repeated  atomic  labels,  that  in  the  third  term  over  four  sets,  and 
that  in  the  fourth  term  over  six  sets.  The  energy  Mr  in  Eq.  (92) 
differs  from  the  strain  energy  per  reference  volume  1If0  used  in 
Sec.  3  by  the  constant  since  continuum  energy 

(31)  vanishes  in  the  reference  state.  Introducing  the  matrix  nota¬ 
tion 


H' 


C ijkl ) : 


:  ffc&i 


d1® 


rrcibc  _  aid  she  c cf_ 

)  ~  d  d  d  ,  d  ,  e 

dr(ij)dr(kl)dP(mn) 


(96) 


A  two-point  elasticity  tensor  (i.e.,  tangent  modulus)  for  homoge¬ 
neous  deformations  can  be  found  as 


dP“A  ^  1  y  <?4>  drU 

dFbB  dFaAdFhB  2H0  i_tj  drc{ij)drdkl)dFaAdFhB 
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Defining  the  finite  strain  measure  P-AB-iP^A  dabFb.B  ~  &Ab)  1 2  and 
noting  that 
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the  average  second-order  elastic  constants  in  the  reference  state 
are  defined  from  atomic  quantities  as 


and  noting  that  dq^l  drbky  =  8b,  the  average  spatial  stress  of  Eqs. 
(90)  and  (93)  becomes 


i*j 


,n  E  r(ij)Hiifklmn)Scd9(kt)Sef<l{mn)  +  ■  ■  ■  (97) 


4  n  tr. 

i*J 
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Since  r^;>=-r the  atomic  stiffness  matrix  in  the  first  of  Eq. 
(96)  exhibits  the  natural  symmetries 


Tiiib  _  Tjab  _  rjba  _  Tjiib  _  rjab  _  rjba  /qq\ 

M(ijkl)  ~  tl{jilk)  ~  n{klij)  ~  ~  M(Jik!)  ~  ~  M(ijlk)  ~  ~  M(klji) 

Similar  symmetries  can  be  deduced  for  H°bjCk]mny  Under  a  homo¬ 
geneous  deformation  q(i7)  =  (F-l)R{,;), 


^=^-E<;>44«><> 

k±l 

+  T7T  E  Htlklmn)Sbd^kl)Sceq\mn)Rtij)  +  •  • •  (") 

k*l 

m+n 

In  the  harmonic  approximation,  products  of  order  higher  than  2  in 
atomic  displacements  q  are  dropped  from  Eq.  (94),  leading  to 


:  <ta,fck^\=,  =  E-E  H(Uk,A,R?iAR(kl) 


"  - 1(ijkl)°.aI^(ij)°.bIS-(kl ) 
0  i+j 
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(103) 


and  clearly  depend  only  on  the  harmonic  part  of  the  potential. 
Similarly,  for  third-order  coefficients 
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Differentiating  Eq.  (102)  with  respect  to  E  gives 
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Third-order  elastic  constants  at  the  reference  state  are  then 
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Relation  (106)  indicates  that  third-order  elastic  constants  depend 
in  part  on  anharmonic  terms  in  the  potential  [48,66].  Derivations 
(99),  (101),  (103),  (104),  and  (106)  rely  on  uniformity  of  the  ato¬ 
mistic  deformation,  such  that  the  relationship  between  current  po¬ 
sition  vector  r  and  reference  position  vector  R  is  linear  in  the 
mapping  F.  In  this  atomistic  context,  F  is  regarded  as  a  linear 
transformation  for  the  primitive  Bravais  lattice  vectors  (e.g.,  akin 
to  Fl  in  Eq.  (12)),  and  though  it  acts  uniformly  over  atoms  in  an 
individual  element  of  reference  volume  fl0,  F  or  FL  need  not  be 
the  gradient  of  any  macroscopic  vector  field  spanning  neighboring 
volume  elements  [51].  In  contrast,  determination  of  deformation 
gradient  measures  for  individual  atoms  of  a  body  undergoing  het¬ 
erogeneous  deformation  is  more  involved  [67],  Expansion  (94), 
while  generic  in  the  sense  that  many  types  of  interactions  (e.g., 
pairwise  and  multibody,  and  central  and  noncentral  forces)  are 
admitted,  may  be  cumbersome  for  computation  of  elastic  con¬ 
stants  for  specific  potentials  not  expressed  explicitly  in  terms  of 
interatomic  separation  vectors.  For  example,  expressions  for 
second-  and  third-order  elastic  constants  for  cubic  crystals  ob¬ 
tained  from  the  embedded  atom  method  are  available  [68]  in 
terms  of  derivatives  with  respect  to  scalar  interatomic  distances. 
However,  these  expressions  [68]  for  second-  and  third-order  elas¬ 
tic  constants  analogously  contain  derivatives  of  up  to  orders  2  and 
3,  respectively,  of  the  potential  functions. 


,-L  y 
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4.4  Net  Residual  Deformation  in  a  Self-Equilibrated 
Lattice.  Now  consider  self-equilibrated  configuration  B  of  Fig.  1 . 
Slip  may  have  occurred  in  achieving  this  configuration  from  the 
reference  state,  but  atoms  occupy  perfect  lattice  sites  apart  from 
the  effects  of  any  defects  remaining  within  the  element  that  either 
(i)  are  also  present  in  the  reference  configuration  B0  or  (ii)  are 
generated  during  the  course  of  plastic  deformation.  Not  consid¬ 
ered  here  is  the  (former)  class  of  defects  present  in  B0.  Although 
positions  of  specific  atoms  differ  in  configurations  B0  and  6  as  a 
result  of  translations  by  Burgers  vectors  of  dislocations  that  have 
passed  through  the  volume  element,  coordinates  R^  can  still  be 
used  to  identify  positions  of  atoms  occupying  perfect  lattice  sites 
in  B  [51].  The  volume  element  may  either  be  treated  as  an  isolated 
ensemble  of  atoms  with  free  boundaries  or  as  part  of  an  infinite 
medium  with  a  periodic  defect  distribution.  The  former  case  is 
more  consistent  with  the  fundamental  definition  of  the  isolated 
unloaded  configuration  [26]  given  in  Sec.  2,  while  the  latter 
case — which  implicitly  includes  the  effects  of  image  forces  of 
defects  in  neighboring  volume  elements — may  be  more  practical 
from  the  standpoint  of  lattice  statics  calculations  with  periodic 
boundary  conditions  [69].  Previous  definitions  for  achieving 
slipped  configuration  B  discussed  in  Sec.  2.1  still  apply  here,  but 


dislocation  fluxes  contributing  to  the  plastic  deformation  are  now 
interpreted  in  terms  of  velocities  of  atoms  comprising  each  line 
defect  rather  than  as  continuum  quantities.  External  forces  vanish 
by  definition  in  B ,  as  do  dynamics  (i.e.,  atomic  velocities  and 
accelerations),  so  lattice  statics  relation  (86)  applies  for  a  self- 
equilibrated  lattice.  However,  because  of  the  presence  of  defects, 

the  stable  total  energy  of  the  ensemble  C  of  atoms  in  B  is  a  local 
minimum  [47,69]  as  opposed  to  the  global  minimum  correspond¬ 
ing  to  perfect  lattice  in  B0.  This  implies  that  some  atoms  occupy 
metastable  positions;  some  interatomic  (i.e.,  internal)  forces  do 
not  vanish  within  the  element  in  B  as  they  do  in  B0  according  to 
Eq.  (95);  and  the  potential  energy  in  B  exceeds  <J>0.  Existence  of  a 
metastable  state  of  local  minimum  energy  results  from  a  lack  of 
strict  convexity  of  the  total  potential  energy. 

The  position  of  atom  i  in  configuration  B  is  denoted  by 

Because  B  is  self-equilibrated,  average  static  atomic  stress  mea¬ 
sures  vanished  by  definition: 


,  1  ^  ,  r)(J> 
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(107) 


with  O,=ja0  the  system  volume  in  B,  <1>  the  potential  energy  of 
Eq.  (87),  and  Because  atomic  coordinates  are  not 

mapped  homogeneously  from  their  positions  in  B0  to  their  posi¬ 
tions  in  B,  the  average  first  Piola-Kirchhoff  stress  measure  in  the 
second  of  Eq.  (107)  does  not  follow  from  a  linear  transformation 
such  as  used  in  the  chain  rule  in  Eq.  (92);  it  is  instead  introduced 
as  a  fundamental  definition.  Equation  (107)  contains  discrete  ana¬ 
logs  of  continuum  relations  (23)  and  (29).  Here  as  in  Eq.  (94),  the 
potential  energy  is  expanded  in  a  series  about  a  perfect  reference 
configuration: 


<£  =  <£(>+  H^ijkI)Sacq^Sbdq^ 
i±j 
k±l 

+  ^  2  Hiijklmn) 3ad<f{i j) ^be^{kl) $cf<f{mn)  +  ■  •  •  008) 

°  i*j 
k+l 
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with  q^=r(i;)-R^  the  displacement  difference  between  two  at¬ 
oms  in  the  defective  lattice,  and  where  Eqs.  (95)  and  (96)  have 
been  used  to  define  atomic  stiffness  coefficients  in  Eq.  (108).  Sub¬ 
stituting  Eq.  (108)  into  the  first  of  Eq.  (107)  and  multiplying 

through  by  J ,  the  null  average  atomic  stress  relation  becomes 

k±l 

+  T7T"  2  4 <llrnn)\/WScMrnn)  +  ■  ■  009) 

0  &j 
k±l 

m+n 


Just  as  in  Eqs.  (36)-(40),  vanishing  of  the  average  symmetric 
stress  measured  in  the  reference  configuration  leads  by  definition 
to  six  equations  for  six  unknown  components  of  the  symmetric 
deformation  map  F1.  Identifying  harmonic  terms  in  the  atomic 
definition  of  the  average  stress  with  linear  terms  in  the  continuum 
elastic  definition  of  the  average  stress: 
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the  following  definition  emerges: 
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with  the  second  equality  following  directly  from  Eq.  (109).  As¬ 
suming  that  the  reference  stiffness  of  the  perfect  lattice  CABCD  is 
positive  definite, 


F CD  ~  &CD 
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where  the  compliance  SCDAB  satisfies  Eq.  (41)  and  can  be  found 
from  inverting  the  second-order  reference  moduli,  the  latter  either 
computed  using  Eq.  (103)  and  the  chain  rule  or  computed  from 
expressions  available  for  particular  potentials  [68],  The  series  in 
brackets  on  the  right  of  Eq.  (112)  accounts  for  anharmonic  inter¬ 
actions,  and  includes  terms  of  order  2  and  higher  in  relative 
atomic  displacements  in  the  locally  deformed,  yet  self- 
equilibrated,  defective  lattice.  For  harmonic  lattice  statics,  in 
which  case  0  and/or  quadratic  and  higher-order  terms  in 

are  dropped,  the  bracketed  series  degenerates  to  zero  and 
Ficd=Scd,  analogously  to  linear  continuum  description  (38).  The 
definition  for  implied  by  the  first  of  Eq.  (Ill),  while  not 
unique,  is  motivated  by  the  corresponding  continuum  definition 
FcD=^lh(c,D)dV=\r'SCDABCABEFfxiEtF)dV  in  Eqs.  (35)  and 
(40).  Following  Sec.  3.6,  the  first-order  accurate  residual  volume 
change  is  7=n/nos=det(4:£)) ~ F!£-(fc+ 1,  where  continuum 
and  discrete  measures  of  reference  and  intermediate  volumes  are 
related,  respectively,  by  fl0=l/  and  ft  =  V+AV.  Calculation  of  this 
volume  change  does  not  require  identification  of  a  bounding  sur¬ 
face  delineating  the  volume  occupied  by  the  atoms  in  configura¬ 
tion  B,  i.e.,  the  atoms  can  occupy  an  arbitrary  shape  so  long  as 
their  interaction  forces  self-equilibrate.  If  such  a  surface  can  be 
identified  in  a  lattice  statics  calculation,  for  example,  an  ensemble 
of  atoms  in  a  defective  state  contained  within  a  hexahedral  bound¬ 
ing  box,  it  should  be  possible  to  compute  FI  trivially  (i.e.,  FI  is 
then  the  volume  of  the  box)  and  then  compare  the  results  with  the 
predictions  of  Eq.  (112),  so  long  as  the  average  static  virial  stress 
vanishes  over  the  defective  ensemble  of  atoms.  An  apparent  ad¬ 
vantage  of  the  atomistic  approach  over  the  continuum  approach  of 
Sec.  3  is  that  effects  of  defect  cores  are  accurately  incorporated  in 
the  former,  presuming  that  the  atomic  potential  is  sufficiently  ro¬ 
bust  to  address  defect  core  structures.  An  apparent  disadvantage  is 
that  solutions  to  Eq.  (112)  must  be  computed  numerically,  while 
the  continuum  elastic  approach  [8,43,44]  provides  convenient 
analytical  formulas  for  volume  changes  such  as  Eqs.  (61)  and 
(74). 


5  Conclusions 

A  three-term  decomposition  for  finite  multiplicative  elastoplas- 
ticity  of  a  volume  element  of  a  crystalline  solid  has  been  de¬ 
scribed.  Isochoric  plastic  deformation  arises  from  displacement 
discontinuities  such  as  those  resulting  from  cumulative  slip  of 
dislocations.  Lattice  deformation  consists  of  two  parts:  (i)  that 
arising  from  recoverable  straining  due  to  external  loading  and 
rotation  of  the  lattice  vectors  and  (ii)  that  arising  from  residual 
elasticity  present  when  the  volume  element  is  self-equilibrated. 
Local  displacement  gradient  fields  within  the  element  contributing 
to  the  latter  (i.e.,  contributing  to  the  average  residual  elastic  strain) 
are  conjugate  to  local  elastic  stress  fields  associated  with  lattice 
defects,  for  example,  eigenstresses  of  dislocation  and  disclination 
lines  and  loops.  Two  methods  are  used  to  quantify  the  average 
residual  elastic  strain.  The  first  treats  the  body  as  a  nonlinear 
hyperelastic  solid  with  spatially  constant  moduli.  In  this  case,  the 
strain  energy  is  expanded  in  a  series  of  up  to  order  3  in  displace¬ 
ment  gradients.  An  integral  equation  for  the  average  residual  elas¬ 
tic  strain  is  obtained  from  the  equilibrium  conditions  requiring  the 
average  stress  to  vanish.  It  follows  that  this  average  residual  strain 
vanishes  identically  only  when  terms  of  order  3  are  discarded,  i.e., 
when  the  body  is  treated  as  linear  elastic.  The  second  method 
follows  from  treating  the  body  as  a  set  of  discrete  atoms  in  the 
context  of  molecular  or  lattice  statics.  The  potential  energy  is 
expanded  in  a  series  about  a  reference  state,  wherein  interatomic 
forces  and  average  atomic  stresses  (i.e.,  average  static  virial  stress 
components)  vanish.  When  the  series  contains  only  terms  of  up  to 
order  2  in  relative  interatomic  displacements  (i.e.,  harmonic  lattice 
statics),  average  dimensional  changes  of  the  lattice  consisting  of 
an  ensemble  of  atoms  in  a  defective  arrangement  are  negligible 
according  to  the  definitions  used  here.  On  the  other  hand,  when 
anharmonic  terms  are  retained  in  the  series  expansion  of  the  po¬ 
tential  energy,  average  dimensional  changes  do  not  always  vanish 
identically.  Particular  attention  has  been  given  to  average  volume 
changes  resulting  from  strain  energy  of  line  defects  in  an  exter¬ 
nally  unstressed  body,  and  examples  for  dislocation  and  disclina¬ 
tion  lines  and  loops  have  been  tabulated  and  compared  with  ex¬ 
perimental  data  for  copper.  Theoretical  predictions  and 
experimental  measurements  for  volumetric  expansion  resulting 
from  dislocation  and  disclination  lines  or  loops  agree  within  about 
a  factor  of  2.  Because  average  residual  elastic  deformation  results 
from  nonlinear  elastic  effects,  i.e.,  products  of  displacement  gra¬ 
dients  of  orders  2  and  higher  in  the  continuum  treatment,  or  from 
anharmonic  effects,  i.e.,  nonlinear  force-displacement  relations  in 
the  atomistic  treatment,  the  contribution  of  average  residual  elas¬ 
tic  deformation  to  the  total  deformation  gradient  will  generally  be 
small  in  conventional  engineering  applications  wherein  defect 
densities  remain  small  to  moderate,  as  observed  experimentally 
for  volume  changes  in  metallic  crystals  deformed  quasistatically 
to  macroscopic  shear  or  compressive  strains  on  the  order  of  unity. 
However,  for  situations  in  which  defect  densities  are  extremely 
large,  such  as  severe  plastic  deformation  processes  or  shock  load¬ 
ing,  the  theory  suggests  that  the  net  contribution  of  residual  elas¬ 
ticity  from  defects,  especially  volume  changes,  to  the  overall  de¬ 
formation  gradient  could  affect  measured  properties  such  as  the 
apparent  bulk  modulus  and  its  pressure  derivatives. 
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